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9^ ABSTRACT 

<: 

■ We treat transformation properties of external radiation fields in the proper frame 
of a plasma moving with constant speed. The specific spectral energy densities of ex- 
ternal isotropic and accretion-disk radiation fields are derived in the comoving frame of 
relativistic outflows, such as those thought to be found near black-hole jet and gamma- 
ray burst sources. Nonthermal electrons and positrons Compton-scatter this radiation 

CN ' field, and high-energy protons and ions interact with this field through photomeson 

■ and photopair production. We revisit the problem of the Compton-scattered spectrum 
associated with an external accretion-disk radiation field, and clarify a past treatment 

■ by the authors. Simple expressions for energy-loss rates and Thomson-scattered spec- 
tra are given for ambient soft photon fields consisting either of a surrounding external 
isotropic monochromatic radiation field, or of an azimuthally symmetric, geometrically 

\ thin accretion-disk radiation field. A model for blazar emission is presented that dis- 

plays a characteristic spectral and variability behavior due to the presence of a direct 
accretion-disk component. The disk component and distinct flaring behavior can be 
X ' bright enough to be detected from flat spectrum radio quasars with GLAST. Spectral 

states of blazars are characterized by the relative importance of the accretion-disk and 
scattered radiation fields and, in the extended jet, by the accretion disk, inner jet, and 
cosmic microwave background radiation fields. 

Subject headings: blazars — galaxies: jets — gamma ray bursts — radiation processes: 
nonthermal 



1. Introduction 

The discovery of intense medium-energy gamma radiation from over 60 blazar AGNs with the 
EGRET instrument on the Compton Observatory (Hartman et al. 1999) shows that nonthermal 
gamma-ray production is an important dissipation mechanism of jet energy generated by black- 
hole accretion. In the standard model for blazars, nonthermal synchrotron emission is radiated by 
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electrons accelerated via first-order shock Fermi processes to Lorentz factors 7 5 x 10 7 /i/i?(G), 
where B is the comoving magnetic field. These electrons will Compton-scatter all available radiation 
fields, including the internal synchrotron field and the external radiation fields intercepted by the 
jet. The intensities of the external ambient radiation fields, transformed to the comoving jet frame, 
are generally required to model the production spectra of blazars. External photon fields include 
the cosmic microwave background radiation (CMBR) field (e.g., Burbidge et al. (1974); Tavecchio 
et al. (2000); Dermer and Schlickeiser (1993), hereafter DS93), the accretion-disk radiation field 
(Dermer et al. (1992), DS93), the accretion-disk radiation field scattered by surrounding gas and 
dust (Sikora et al. 1994; Dermer, Sturner, and Schlickeiser 1997), infrared emissions from hot 
dust or a molecular torus (Protheroe and Biermann 1997; Blazejowski et al. 2000; Arbeiter, Pohl, 
and Schlickeiser 2002), reflected synchrotron radiation (Ghisellini and Madau 1996; Bottcher and 
Dermer 1998), broad-line region atomic line radiation (Koratkar et al. 1998), etc. 

Gamma rays in leptonic models of broadband blazar emission originate from synchrotron self- 
Compton (SSC) (Maraschi, Ghisellini, & Celotti 1992; Bloom & Marscher 1996; Tavecchio et al. 
1998) or external Compton (EC) (e.g., Melia and Konigl (1989), DS93, Bottcher, Mause, and 
Schlickeiser (1997); see Bottcher (2001) and Sikora &: Madejski (2001) for recent reviews) pro- 
cesses. In hadronic models, secondary photopairs and photomesons are produced when energetic 
protons and ions interact with ambient synchrotron photons (e.g., Mannheim and Biermann (1992); 
Mannheim (1993)) and photons of the external field (Bednarek and Protheroe 1999; Atoyan & Der- 
mer 2001). The spectral energy distribution of a blazar will consist, in addition to nonthermal 
synchrotron and SSC emissions, of secondary radiations that are a consequence of processes involv- 
ing external soft photon fields. These radiations include neutral secondaries from hadrons, and the 
radiations produced when external soft photon fields are scattered by directly accelerated lepton 
primaries and lepton secondaries formed in photopair, photohadronic, and electromagnetic cascade 
processes. The intensities of the various components depend on the properties of the external fields, 
the properties of the relativistic outflows, and the time-dependent spectral injection of ions and 
electrons into the outflow. Given these conditions, the evolving particle and photon distributions 
and emergent radiation spectra can be derived for a specific model. 

A significant modeling effort has been devoted to blazar studies. Where the synchrotron and 
SSC components are the dominant radiation processes, as appears to be the case in some X-ray 
bright BL Lac objects such as Mrk 421 and Mrk 501, observations of correlated X-ray/TeV loop 
patterns found in graphs of spectral index versus intensity are explained through combined particle 
acceleration, injection, and radiative loss effects (e.g., Kirk, Rieger, and Mastichiadis 1998; Li and 
Kusunose 2000). The spectral energy distributions of BL Lac objects are successfully modelled with 
synchrotron and SSC components from broken power-law electron distributions with low energy 
cutoffs, as in Mrk 501 (Mastichiadis and Kirk 1997; Pian et al. 1998). Models for the spectral 
energy distributions of flat spectrum quasars, including external Compton components, have been 
presented for different spectral states by Hartman et al. (2001), Mukherjee et al. (1999), and 
Bottcher (1999), and in papers cited in the first paragraph. 
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More recently, Sikora et al. (2001) have presented a model for blazar variability with time- 
dependent injection into a relativistically moving blob that contains nonthermal electrons which are 
subject to radiative and adiabatic losses. The scattered radiation from an external, quasi-isotropic 
radiation field is treated, and evolving spectral energy distributions and light curves are calculated. 
No direct accretion disk field is however treated, though the emission sites could be within 0.1-1 pc of 
the black hole, where such components can make a significant contribution (Dermer and Schlickeiser 
1994). We show that the addition of the accretion-disk components produce light curves and 
correlated multifrequency behaviors that are quite distinct from the behavior calculated by Sikora 
et al. (2001). These distinct patterns are bright enought to be well detected with the upcoming 
Gamma ray Large Area Space Telescope (GLAST) mission, 1 thus constraining the location of the 
shocks in gamma-ray blazars and the dominant spectral components. 

In this paper, we revisit the transformation properties of the external radiation field, and 
clarify a past treatment by the authors (DS93). Here we adopt a formulation in terms of intensity 
rather than, as in DS93, photon number. The intensity formulation seems to offer a simpler method 
for comparing multiwavelength data, though the photon number approach perhaps simplifies the 
formulation of the scattering, attenuation, and secondary production interactions. 

The jet plasma is assumed to be located along the symmetry axis of a geometrically-thin 
accretion disk. The intensity of photons in the jet frame depends on whether the accretion disk 
is optically thick or thin at large disk radii R 3> r, where r is the distance of the jet plasma 
from the black hole. The treatment of DS93 used an optically-thin formulation with a volume- 
energy generation rate appropriate to disk accretion in a steady flow. DS93 assigned a blackbody 
temperature for the mean energy of the disk radiation field at different annuli, whereas a wide 
variety of models are possible. The radiation fields originating from radii R » r, though spectrally 
distinct between the optically-thin and -thick cases, give minor effect on the Thomson energy losses 
and scattered photon spectra. 

The energy density of disk photons that are scattered to 7-ray energies by the jet electrons 
is shown to be dominated by photon production at R ~ r in the near-field (NF) regime, and by 
the total black-hole power approximated as a receding point source in the far-field (FF) regime. 
The accretion-disk photons scattered by gas and dust surrounding the black hole can dominate 
the near- and far-field disk fields at large gravitational radii (r g = 1.5 x 10 14 Mg cm), depending 
on the relative transition radii where the NF, FF, and scattered external radiation fields dominate 
(Appendix A). 

The transformed energy spectrum of the thermal accretion disk in the comoving frame displays 
a low-energy steepening for photons which originate primarily from radii R S> r. This is an 
important effect in photohadronic calculations of very high energy secondary production, and is 
also important in 7-7 transparency calculations (Atoyan & Dermer 2002). Other photon sources can 



1 see http://glast.gsfc.nasa.gov 
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make an additional contribution in this range. Because the bulk of the black-hole accretion power 
is generated at small disk radii, the use of an optically-thin or optically-thick formulation at R S> r 
has negligible effect on the total scattered power for power-law electron injection with injection 
index > 2. This article presents basic equations to calculate transformed external radiation fields 
from geometrically-thin accretion disks in the limiting optically-thin and optically-thick regimes. 

Expressions for Thomson-scattered spectra from nonthermal electrons with arbitrary energy 
and isotropic comoving pitch-angle distributions are derived for external isotropic and direct accretion- 
disk radiation fields. The electron equation-of- motion is solved in the regime where the total syn- 
chrotron and external Thomson losses from the disk and quasi-isotropic external radiation fields 
dominate SSC losses. The evolving electron distribution is obtained through a continuity equation 
treatment with radiative losses on stationary electron injection over a range of radii, but here ne- 
glecting adiabatic losses. This simplified blazar model shows a dramatic change of spectral behavior 
if nonthermal particle acceleration continuously operates within £rf 10 2 -10 4 r 9 from the black-hole 
engine. In the model studied here, a pivoting of spectral shape is expected at a few GeV on time 
scales of days, with 7-ray spectral activity correlating with monotonic increases at IR and X-ray 
energies. Delayed radio activity is predicted as the emergent energized plasma blob expands and 
becomes optically thin to synchrotron self-absorption on the parsec scale, though this behavior is 
not modeled here. Campaigns employing GLAST accompanied by broadband monitoring at other 
frequencies can record multifrequency spectral flaring data and test the model. 

Preliminaries of radiation theory and invariants are given in Sections 2 and 3, respectively. 
Section 4 treats the transformation of external radiation fields, specializing to the azimuthally 
symmetric and extreme-relativistic (r> 1) cases. The comoving specific spectral energy density 
for an optically-thick Shakura and Sunyaev (1973) accretion disk is derived in Section 5. Accretion- 
disk modeling is briefly mentioned (§5.1) before proceeding to derive the comoving energy spectrum 
of the disk field (§5.2). The same problem for an optically-thin disk is treated in Section 6, and 
comparison with the approach of DS93 is made in Section 7. The properties of the specific spectral 
energy density of the accretion-disk radiation field, transformed to the comoving fluid frame of the 
jet, are examined in Section 8. 

Section 9 presents analytic forms for the energy-loss rates and spectra from various radia- 
tion processes, including synchrotron emission (§9.1), radiation from Thomson-scattered external 
isotropic monochromatic (§9.2) and accretion-disk fields in the near-field (§9.3) and far-field (§9.4) 
limits, and SSC radiation (§9.5). A simple model for blazar variability is presented in Section 
10 in the limit that synchrotron and Thomson-scattered external radiation fields dominate the 
electron-energy loss rates. Model results for stationary particle acceleration and injection tak- 
ing place within ~ 10 2 -10 4 r 9 from the black hole reveal a characteristic spectral behavior that 
joint radio/submillimeter/IR/X-ray/7-ray campaigns can test. Detectability of blazar flares with 
GLAST is treated in Section 11, where multiwavelength light curves of a model blazar flare are 
presented. 
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The Appendices provide more detail on the derivations of the results in Section 9. In particular, 
Appendices A. 6 and A. 7 show that different spectral states of blazars and extended jets are defined 
by the simple relations in equations (A13) and (A15) which govern the relative imporatance of the 
accretion-disk, scattered, and CMBR fields. An important deficiency in the model is the treatment 
of Klein-Nishina effects in Compton scattering, which has been considered in the context of blazar 
physics by Bottcher, Mause, and Schlickeiser (1997), Georganopoulos, Kirk, & Mastichiadis (2001) 
and Dermer and Atoyan (2002). An analytic approach to this physics is given in Appendix D, 
including comparison with the paper by Georganopoulos, Kirk, & Mastichiadis (2001). 



2. Preliminaries 

We consider external radiation fields from optically-thick and optically-thin accretion disks, 
following preliminaries on radiation physics (Novikov and Thorne 1973; Rybicki and Lightman 
1979). Our notation is such that a quantity X t has the full units of A, whereas differential quantities 
are denoted by X(xi,x 2 ;y) = dX(y)/dxidx2, where y is a parameter. 

The intensity 

= dAdtdedSl (1) 
is defined such that I e (Sl)dAdtdSl is the infinitesimal energy d£ carried during infinitesimal time dt 
by photons with dimensionless energy e = hv/mg.c 2 between e and e + de that passes through area 
element dA oriented normal to the direction of rays of light lying within the solid angle element 
dSl about the direction SI. Except for polarization, I £ (Sl;x,t) gives a complete description of the 
radiation field at location x and time t. 

The quantity dF e = [iI e (Sl)dSl is the differential spectral energy flux (units of ergs cm" 2 s -1 
e _1 ) that passes through an area element oriented at an angle /x = cos = n ■ SI to the rays within 
solid angle element dSl about the direction SI. The net flux F £ (n) = § dSl nI £ (SX) depends on the 
orientation of the area element. Radiation fields that are azimuthally symmetric about some axis 
and that can be expanded in even powers of jjl have zero net flux, that is, F £ = 0. Obviously, the 
net flux of an isotropic radiation field is zero. 

Because the momentum of a photon is hv/c, the component of momentum flux along h is 
given by p e (dynes cm _2 e _1 ) = c~ l § dSl fi 2 I e (Sl). One factor of /u is from the projection of the area 
element, and the second is from the component of momentum along h. Although I e (Sl) > for all 
p e can be < if the net momentum flux is directed opposite to h. The total radiation pressure 
is p = c -i J o °° de § dn n 2 i £ (n). 

The vF v energy flux (ergs cm~ 2 s _1 ) at frequency v = m e c 2 e/h is / e (n) = e F e = e § d^l pI e (Q). 
The bolometric energy flux 3>(n) = f^°dee~ 1 f e (h) = de§ d£l fiI £ (Q,). The specific spectral 
energy density u(e, SY) = dS/dVdedQ, and the spectral energy density u(e) = <fd£lu(e, SI). Note 
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that dS = u(e, Q)dVdedQ = u(e, £l)cdtdAded£l = I e (Q)dAdtdedQ, so that 

u(e,n) = c - 1 i e (n) . (2) 

The mean intensity I e = (4vr)~ 1 § dtl I e (tl) = cu(e)/(4vr). Thus u{e) = 4irl e /c. The total 
radiation energy density u = J °° de u(e) = Attc^ 1 f °° de I e . The radiation pressure for an isotropic 
radiation field with (mean) intensity I l e so is p = c" 1 / °° de § dfl fi 2 !™ = (4ir/3c) J °° de I l e so . Be- 
cause u = (4ir/c) / °° de I* so , we have 

p = u/3 (3) 

for an isotropic radiation field. 

From conservation of energy, we can say that dE\ = I e ^dAidtd^l\dei = I 6i2 dA 2 dtdQ, 2 de 2 by 
following a bundle of rays from one area element to the next. The relationships between area 
elements separated by distance d, and the solid angles subtended by the bundle of rays at the two 
locations are dfl\ = dA 2 /d 2 and dfl 2 = dAi/d 2 . With e\ = e 2 = e, de\ = de 2 = de, we obtain 

4,1 = 4,2, or dljds = 0, (4) 

provided no absorption occurs as the light rays propagate. The equation of radiative transfer is 

^ = ~ K eh + i(e, 0.) , or ^ = -I e + S e , (5) 



ds dr. 



c 



where we define the photon-energy dependent spectral optical depth r e = n e ds in terms of the 
e-dependent spectral absorption coefficient K e (cm -1 ) and the emissivity j(e, SY) = d£ /dVdtdedQ. 
The source function is 

S e = J ^. (6) 

The energy flux from a uniform brightness sphere of radius R is an instructive example. In 
this case, / e (fi) = B e , where B e is independent of Q by definition. Further defining J °° deB t = B, 
we find that the total energy flux at the telescope aperture located a distance d from the source is 
<3?(n) = J °° de§ dft [iB e = 2ttB f^/ 1 _ R 2/ r 2 ^ dfi = nBR 2 jr 2 . Thus the energy flux at the surface 
of a uniform brightness sphere is irB. The energy flux at the surface of a blackbody radiator is 
usbT 4 , where the Stefan-Boltzmann constant ctsb = 5.67 x 10~ 5 ergs s _1 cm~ 2 degrees" 4 . 



3. Invariant Quantities 

The elementary invariants are the invariant 4-volume d 3 xdt = dVdt, the quantity d 3 p/E, and 
the invariant phase volume dV = d 3 xd 3 p. Here E is photon or particle energy, p is momentum, 
and d 3 p = p 2 dpd£l. The invariance of d 3 xdt and d 3 p/E can be demonstrated by calculating the 
Jacobian for transformations consisting of rotations and boosts. The invariance of dV follows by 
noting that dt/E is the ratio of parallel 4- vectors (Blumenthal and Gould 1970). 
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Because the number N of particles or photons is invariant, 

— - J- dN - 1 dN 1_ dE _ _ 3 u(e,Q) 

~dV ~p> dVdpdfl ~ (m e c) 3 e 2 dVdedf} ~ mjc 5 e 3 dVdedtt ~ 6 mfc 5 ' ^ 

where the latter three expressions apply to photons (E = mc 2 e). Hence, e~ 3 u(e,Q) and e~ 3 / e (^) 
are invariants. 

The function 

dN = 1 rf£ = 1 ., , , . 

d 3 xdtd 3 p {m e cfe 2 dVdtdedn (m e c) 3 e 2 ^ J 1 j 

is invariant, where the last two expressions apply to photons, implying that e~ 2 j(e,$7) is invariant. 
This and the invariance of e~^I e (Q) show from equation (6) that en e is invariant. 



4. Transformation of External Radiation Fields 

The photon energy e and angle 9 = cos -1 fi in the stationary frame of the accretion disk are 
given in terms of the comoving photon energy e' and direction cosine fi' through the relations 
e = re'(l + /?r/-0> A* = (/•*' + A^r)/(1 + Prv'), and (f> = <j)' . From the invariance of u(e, ri)/e 3 , we have 

For simplicity, we assume azimuthal symmetry about the axis of the boost to the comoving jet 
frame, moving with bulk Lorentz factor T with respect to the stationary black hole system. For 
an isotropic external photon field, u(e,fi) = \u(e). Hence u'(e',n') = n(e)/[2r 3 (l + (3r^') 3 ] in this 
case. 

For an external isotropic monochromatic radiation field, u(e) = u ext 5(e — e). Thus the comoving 
energy density 

u' ext = J™ de> J 1 d»>u'(e', „') = g§ J 1 V (1 + /? r //r 4 = u ext r 2 (l + (10) 

(Dermer and Schlickeiser 1994). Following the integration over energy, one sees that the angle- 
dependent specific energy density u'(n') = u ex t/[2T A (\ + (3yh') a }- The four powers of the Doppler 
factor 5 = T(l + flrn') = [T(l — /^rA 4 )] -1 are due to two powers of S from the solid angle transforma- 
tion, one from the energy transformation, and one from increased density due to length contraction. 
For r> 1, the specific energy density ranges in value from ~ 8T A u ex t at jj! = — 1 to ~ u ex t/(32T 4 ) 
at fi' = 1. The function plummets in value when /j,' > — (3r- Multiplying u'(fj! = —fir) = T A u ex t/2 
by the characteristic solid angle element 5Q' ~ tt6' 2 ~ tt/F 2 gives u' ext = (it /2)T 2 u ext , which approx- 
imately recovers equation (10) when T > 1. The function u'(e', //) = 4T 2 u ext 5(e' -4Te/3)5(fi' + l)/3 
provides a useful approximation in the limit f » 1, noting that u' ext /n' ext = (1 + f3^/3)Tem e c 2 . 
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For isotropic external photons, the hardest and most intense radiation is directed opposite to the 
direction of motion of the jet plasma in the comoving frame. 

Now consider the transformation of an isotropic, azimuthally symmetric power-law radiation 
field described by the function 

u(e,fi) = k a u ext e~ a , for e e < e < e u . (11) 
The normalization u ex t = Jq° de J^ 1 dfiu(e, n) gives 



k a 



9 , i-i "I., , for general a 
^[2\n{e u /e t )}-\ for a = 1 . 



Equations (9) and (11) imply 



(12) 



^ ^ = v\i+t7? +a ' for et < re ' (i + M < eu ' (13) 

which is easily shown to be normalized to u' ext given by equation (10). 



5. Specific Spectral Energy Density 

The specific spectral energy density u(e,Q;x, t) = c _1 /(e, ft) from equation (2). The external 
field at arbitrary locations can be obtained from the constancy of specific intensity, equation (4), 
if intervening absorption is negligible. This is the starting point for the calculations of external 
radiation fields, from which transformed fields can be obtained through equation (9). 

We specialize to the calculation of the external radiation field emitted by steady, azimuthally 
symmetric emitting regions, as might correspond to thin accretion disks around supermassive black 
holes. We furthermore consider a geometry where the emitting surface of the disk is located on the 
symmetry plane, and assume that the disk thickness h(R) at radius R is <C R- Thus 

= (l+R 2 /r 2 )- 1/2 , R = r(fi- 2 - 1) 1/2 • (14) 



Vr 2 + R 2 



In the cool, optically-thick blackbody solution of Shakura and Sunyaev (1973), the disk emis- 
sion is approximated by a surface radiating at the blackbody temperature associated with the 
local energy dissipation rate per unit surface area, which is derived from considerations of viscous 
dissipation of the gravitational potential energy of the accreting material (Shapiro & Teukolsky 
1983). One finds that the integrated emission spectrum measured far from the black hole varies 
oc e 1 / 3 up to a maximum photon energy associated with the innermost stable orbit of the accretion 
disk. The optically-thick solution is unstable in the inner region near a black hole due to secular 
density/cooling instabilities (Lightman and Eardley 1974) in some regimes of the Eddington ratio 

^Edd = t , (15) 

-^Edd 
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where rj is the efficiency to transform accreted matter to escaping radiant energy. The Eddington 
luminosity LEdd = 1.26 x 10 47 Mg ergs s _1 , where the mass of the central supermassive black hole 
is M = W 9 MgM Q and the black hole is accreting mass at the rate m (gm s -1 ). 



The intensity of a blackbody is 
2»(0) = 2c(m e c 2 ) e 3 2c(m e c 2 ) 



A 3 [exp(e/G) - 1] A 3 C 



' ©e 2 , £<9 Rayleigh-Jeans regime 

, e 3 exp(— e/0) , e>9 Wien regime 



(16) 

where = kBT/m e c 2 and Ac = h/m e c = 2.426x 1CT 10 cm is the electron Compton wavelength. For 
steady flows where the energy is derived from the viscous dissipation of the gravitational potential 
energy of the accreting matter, the radiant surface-energy flux 

dE 3GMm 

lAlt = "W < 17 > 

(Shakura and Sunyaev 1973), where 

<p(R) = [l-0 i (R i /R) 1/2 }, (18) 

Pi = 1, and Ri = 6GM/c 2 for a Schwarzschild metric. Integrating equation (17) over a two-sided 
disk gives 77 = 1/12. Assuming that the disk is as an optically-thick blackbody radiator, the effective 
temperature of the disk can be determined by equating equation (17) with the surface energy flux 
a SB T 4 (i?). 



5.1. Disk Models 

To calculate the scattered jet radiation spectrum from the disk-jet system, it is essential to 
properly characterize the accretion-disk geometry and emissivity. Observations suggest that the 
radius separating a geometrically-thin outer accretion disk from an optically-thin hot inner cloud or 
disk increases with decreasing ^Edd in the range ^Edd i$0.1. (see Liang (1998) for a review of accretion 
disks around galactic black hole sources). Within an advection-dominated accretion disk (ADAF) 
scenario (e.g., Esin, McClintock, & Narayan (1997); Di Matteo et al. (2000); Quataert & Gruzinov 
(2000)), the radiant luminosity from an accreting system declines markedly when ^Edd C0.1 due 
to the advection of photons into the black hole and the convection of angular momentum and mass 
outward due to convective instabilities in the advection-dominated flows. Thus rj = r/(m) in eq. 
(15), and the radiant luminosity follows a steeper dependence than L rac i ^5 -^Edd when ^Edd *C 0.1. 
The radio-loud branch of accreting black holes seems to be found in the convectively unstable, low 
Eddington luminosity regime. 

In this paper, we consider only the simplest flat-disk spectrum, which is assumed to be well- 
approximated by the Shakura/Sunyaev disk spectrum. Optically-thin emission from the inner 
disk may also be present and can be treated according the formulation presented here, but is 



-10- 



not considered here. The ejection of relativistic jet plasma probably occurs in the low-luminosity 
(^Edd ^5 0-1) regime, and evolutionary considerations (Bottcher and Dermer 2002) are consistent 
with this inference. 

The flux from this disk model is given by 

(Shapiro & Teukolsky 1983). Here we use a monochromatic approximation for the mean photon 
energy, with T(R)= [3GMrhip(R) /SirR 3 ^^ 1 ^ . It is simple to calculate the transformed radiation 
field in the comoving frame of relativistic plasma along the symmetry axis of the disk with this 
expression. 

The outer accretion disk could be optically thin to Thomson scattering, as we show in Section 
6. When this happens, it is no longer acceptable to adopt an optically-thick scenario. Most of 
the energy is dissipated in the central regions, because AE oc AR/R 2 . Reprocessed radiations can 
dominate viscous radiations in the outer disk, especially for flared disks. The intensity of the outer 
disk can then be dominated by reprocessed UV and X/7 radiation, or emission from a surrounding 
torus. In any case, the emission of the outer disk is unlikely to be given precisely by equation (19), 
though it provides a useful functional form for further study. 



5.2. Integrated Emission Spectrum for Blackbody Disk Model 

The spectral energy density u(e) = c _1 <f cZQ lf°(Q; R) along the jet symmetry axis is evaluated 
for the Shakura/Sunyaev optically thick blackbody disk model from equation (19). We define 
K = 2.7/c B (3GMm/8^s B ) 1/4 /m e c 2 . Hence 

„ ss(£) _ f~ 4, m S[e _ K^ { R)R-^ (20) 

Eq. (20) can be solved analytically in the approximation ip(R) = 1, which becomes accurate in the 
limit e <C K/R 3 ^ 4 , giving 

3GMmr f 00 5(e - KR- 3 / 4 ) _ GMmr (e/K) 1 / 3 
USs{ - e > ~ 8vrc Jr. R 2 (r 2 + i? 2 ) 3 / 2 ~ 2ttKc [r 2 + (K/e) 8 / 3 ] 3 / 2 ( " ' 

GMm \ ' for ^ « ^ « R ^ 

[ r(f ) 13 /3 j for |« r -3/4 

Thus we see that the rapid decline in the energy radiated in the blackbody disk, compounded by 
the small solid angle subtended by the area of the disk when R^> r, leads to a spectral steepening 
in the integrated energy density of the external radiation field. 
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Scaled quantities are adopted in order to get a better idea of the photon energy of the external 
radiation field measured along the disk axis for the blackbody disk model. The mass accretion 
rate m = ^Edd^Edd/(^?c 2 ) = 1.4 x 10 26 £Edd-^9/?/ g m s_1 ) using equation (15). We write r = fr g 
and R = Rr g , where f is the jet height and R is the disk radius in units of gravitational radii 
r g = GM/c 2 = 1.48 x 10 14 M 9 cm. The characteristic photon energy emitted from the disk at R 
gravitational radii from the nucleus is 

meC ^ * 77 (^d )V4^-3/4 eV (22) 

Replacing the term R in equation (22) with the dimensionless jet height f gives the photon energy 
at which the integrated spectral energy density «ss(e) displays a break from the e 1//3 spectrum to 
the e 13//3 spectrum (eq. [21]). For a 10 9 M Q black hole accreting at the Eddington limit at efncency 
rj = 0.1, we therefore see that the mean photon energy radiated at 10, 10 2 , 10 3 and 10 4 gravitational 
radii from the central source is 26, 4.5, 0.81 and 0.14 eV, respectively. The photon energy derived 
at 10r 9 is, however, not that accurate because of the approximation <p — > 1. Other effects such as 
gravitational redshifting also become important when R < 10. In any case, a modified blackbody 
or optically-thin disk model rather than a blackbody model may hold in the inner disk region. 
Evidence for a blackbody Shakura/Sunyaev disk around supermassive black holes is provided by 
the intense optical/UV "big blue bump" radiation, as observed for example in the UV spectrum of 
3C 273 (Lichti et al. 1995; Kriss et al. 1999). 



6. Transformation Properties of the Optically Thin Disk Radiation Field 

The continuity of mass in a steady flow consisting of ionized hydrogen implies a mass accretion 
rate m = 2tt Rh(R)v rat j i (R)m p n(R) , where h(R) is the full disk thickness at radius R, v ra d(R) = 
cflrad is the radial flow speed, and n(R) is the proton density at R. The vertical Thomson scattering 
depth through the disk is 

t t = a T h(R)n(R) = ^ 2.0 - ^ Edd . , (23) 

V ' y ' 27rRc(3 rad m p Rprad(R)v 

If the middle-outer disk rotates in Keplerian motion, than mv 2 ~ GMm/R implies that (3$ ~ R~ l l 2 , 
where v$ = f3gc is the azimuthal speed of the accretion flow. Letting (3 ra d = k(3$ (k < 1) implies 
tt — 2i-^^/krjR 1 / 2 . If the radial flow speed is sufficiently rapid, that is, if R > (2£Edd/^ r /) 2 ) then 
the accretion disk is optically thin to Thomson scattering. This condition requires rj <C iEdd , which 
is compatible with ADAF models (§5.1). 

Within the Newtonian approximation for a thin accretion-disk geometry (Shapiro & Teukolsky 
1983), the surface energy flux is given by equation (17). If the emitting region is optically thin to 
Thomson scattering, then dE/dVdtdVL = [2Trh(R)]~ 1 (dE / dAdt) and the emissivity is given by 

^■•z^T^Ktm^-^* (24> 
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where a monochromatic approximation to the emission spectrum is made. The equation of radiative 
transfer (5) for an optically-thin region gives, noting that As = h(R)//j, for a thin disk, 

R) = T^Wil v{me " m] ■ (25) 



7. Comparison with the Approach of Dermer and Schlickeiser (1993) 

The direct disk radiation field provides the most intense radiation field to be scattered to 
gamma-ray energies when the relativistic plasma ejecta is within some 10 3 -10 4 r 9 from the cen- 
tral source (Dermer and Schlickeiser 1994), and even farther if the scattered radiation is weak. 
To calculate the external radiation field intercepted by the ejecta, DS93 considered the emission 
function 

n(e,n,R) = ^ 5 -^,x= r -, (26) 
4ttx z c 2tt h 

which is related to the specific spectral energy density according to 

roo 

u(e,tt) = m e c 2 dRen(e,n,R). (27) 

JRr 

Equation (26) allows a large range of radially-dependent emissivity functions to be tested, though 
the formulation assumes optically-thin emissivity. 

Note the relations 

em e c 2 N(e,R) = E(e, R) = 2 x 2nR x x 5[e - e(R)} = <p(R) 5[e - e(R)}. (28) 

Therefore 

3GMm [°° Sfr-jiMR) S[e - i(R)\ 

UDs(e ' 0) = 4 dR ^ • (29) 

Transforming the 5- function in fx to a 5- function in R using equation (14) gives 

UBs(e,n) = <p(R)S[e - e(R)] , R = ry^^T, (30) 

equivalent to equation (25). By integrating equation (30) over in the approximation that <p(R) ~ 
1, equation (21) is recovered exactly when r~ 3 / 4 <C e/K <C R] 3 ^ 4 . When e/K <C r -3 / 4 , the 
exponent 13/3 for the optically-thick result becomes 3 in the optically-thin formulation (compare 
eq.[21]). 

The approach of DS93 considers the external radiation field from an optically-thin accretion 
disk that emits photons at a temperature or energy corresponding to the blackbody value. Other 
optically-thin disk models can be developed, including one- and two-temperature disk models, or 
hybrid models, so we treated a specific case in our paper. This description is valid to determine the 
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transformation properties of the disk when treated as an optically-thin radiator with a blackbody 
disk spectrum, and differs only marginally from a blackbody except when scattering photons with 
energies e < A>~ 3 / 4 . The spectral effect makes little difference for Compton scattering rates and 
spectra, but can be important in processes with thresholds, e.g., photopion mechanism and 7-7 
absorption (Atoyan & Dermer 2002). 

A more general formulation must consider the transition between the optically-thin and optically- 
thick regimes, geometrically-thick accretion disks and more general disk models, and emergent jet 
physics. 



8. Transformed Energy Density of the Accretion Disk Radiation Field 

As an illustrative example of the results presented in this paper, we consider the optically- 
thick, geometrically-thin accretion-disk radiation field (eq. [19]) in the approximation ip — ► 1, for 
which the specific spectral energy density along the jet axis is 

"I'-v-mm*-*- <31> 

The energy loss rate in the Thomson limit depends, for a distribution of electrons with random 
pitch angle, only on the total comoving frame energy density 



u' 



ext 



/ da'u'(u') = <k dU de' ( — YVe.fi) = 5— / da' n , ,„ rr , 

(32) 



Fig. 1 shows the integrand of the rightmost integral in equation (32). When r> 1, two domi- 
nant components of the differential energy density make up the total energy density: a component 
from the disk radiation field at disk radii R rs r, called the near-field (NF) component; and a far 
field (FF) component coming directly from behind, which dominates the disk contribution at large 
radii (Dermer and Schlickeiser 1993). The accretion disk radiation field, when approximated as a 
point source that illuminates the ejecta blob directly from behind, presents a total comoving energy 
density 

1 3GMrh 2 

Ub ~ P(i + ^)2 (1 " Y l) ■ {66) 

This result can be derived from the relation Ub = 2-ir J^ max dfi f^°de(lf s /c), using equations (19) 
and (18) in the limit r 3> R, combined with the relativistic transformation of the energy den- 
sity (eq.[6] in Dermer and Schlickeiser (1994)). The point-source approximation improves as the 
accretion disk looks more like a point source in the comoving frame, that is, when 

» max ~ , _ ^ _ Kj2r2) ~ 2( r ) -1. (34) 
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The far- field approximation thus holds only when r 3> TR m i n , and dominates the near- field 
component only when r 3> T 4 R m i n (Appendix A. 3). The following relationship connects the disk 
model with the disk radiant luminosity, neglecting advective effects: 

f°° , n r, dE 3GMm 2ft, a t GMm^ . , 

In Appendices Al and A2, analytic properties are derived of the integrand of the right-most term 
in equation (32) in the NF and FF limits, respectively 

9. Analytic Energy-Loss Rates and Spectra 

This provides sufficient theory to derive expressions for radiation spectra in blazar jets. Sim- 
plified analytic expressions for electron-energy loss rates and spectral components are presented in 
the near- field and far- field regimes. For completeness, we also summarize results of our previous 
work for synchrotron and synchrotron self-Compton processes (Dermer and Schlickeiser 1993; Der- 
mer, Sturner, and Schlickeiser 1997). These expressions employ 5-function approximations for the 
emission spectra, and the spectral forms are least accurate near endpoints and spectral breaks. 

9.1. Synchrotron Radiation 

The rate at which a randomly ordered pitch-angle distribution of relativistic nonthermal elec- 
trons lose energy via the synchrotron process in a region with mean comoving magnetic-field inten- 
sity B is given by —m e c 2A f syn , where 

-isyn = \ ^ 7 2 = 1-3 x 1(T 9 5 V s- 1 = W (36) 
3 m e c z 

(Blumenthal and Gould 1970), and ub = B 2 /8ir is the energy density of the magnetic field. 

In the (5-function approximation for the elementary synchrotron emissivity, the vF v synchrotron 
radiation spectrum from a uniform blob, which is assumed to be spherical and to have a randomly 
oriented magnetic field in the comoving frame, is 

fr ^ ( ^ NM ,^^±. (37) 

In this expression, emission properties are integrated over variations on the comoving size scale 
of the plasma blob. On the corresponding observer time scale t var = (1 + z)rj,/cS, the blob upon 
reaching location x is assumed to host an instantaneous electron energy spectrum N^j). Here 
N e (j; x)d / y is the differential number of electrons with comoving Lorentz factors between 7 and 
7+ d7- N e (j;x) is evaluated by solving a continuity or diffusion equation. Better descriptions of 
the system require integrations over the emitting volumes that arise from light-travel time effects 
(Chiaberge and Ghisellini 1999; Chiang and Dermer 1999; Granot et al. 1999). 
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9.2. Thomson-Scattered External Isotropic Monochromatic Radiation Field 

Consider the Thomson energy-loss rate of electrons that are randomly distributed throughout 
a plasma blob that passes through a uniform external isotropic monochromatic radiation field with 
mean photon energy e* = 10 _4 e_4. The energy-loss rate due to Compton-scattered CMBR in the 
Thomson limit is given through 

-j T = t ^Z_ 7 Vr 2 (l + ^ 33 x 10 -8^[^_] 2 r 2 (1 + /? 2 /3) g -l ^ k 2 (3g) 

3 m e c z cm^ 

where is the energy density of the external isotropic radiation field photons which are scattered 
in the Thomson limit (Dermer and Schlickeiser 1993, 1994). Equation (38) holds for 7 <C 7kn, 
where 7kn is the Lorentz factor where Klein-Nishina (KN) effects become important, given through 
47KNre* = 1. The transition from Thomson scattering to KN scattering takes place at observed 
photon energies e K N = <V[ 16re *(! + z)\ = 300(<5/T)[e_ 4 (l + z)]~ l MeV (see Bottcher, Mause, and 
Schlickeiser (1997); Georganopoulos, Kirk, & Mastichiadis (2001), and Dermer and Atoyan (2002) 
for treatments of KN effects on external Compton scattering) . 

The vF v spectrum radiated by relativistic electrons which scatter photons from an external 
isotropic monochromatic radiation field in the Thomson limit is 




(Dermer, Sturner, and Schlickeiser 1997). The mean ambient CMB photon energy e* = e(l + z), 
e = 2.1k^,TcMBlfne ( ? = 1-24 x 10~ 9 , where Tqmb = 2.72 K is the present temperature of the 
CMB, and u* = 4 x 10~ 13 (1 + z) 4 ergs s _1 is the local ambient CMBR energy density. For blazars, 
the ambient photon field might be reprocessed UV accretion disk radiation, with energy density 
depending on principle resonance line transitions in the broad line region or illuminated disk and 
torus. A simple prescription for the scattered radiation field is to let 

u * = Z^IT = 3 - 5 x 1Cr4 (777^r) > where L d = £ E dd^Edd = r]riic 2 ergs -1 , (40) 

4irrj c c rj c (pc) 0.01 

t sc is an effective scattering depth of the surrounding medium including the broad line region, 
and r sc is the characteristic size of the surrounding scattering medium (Dermer and Schlickeiser 
1994). Effects of spatially varying external radiation energy densities from scattering media with 
power-law density gradients are treated by Blandford & Levinson (1995). 

The distance along the jet axis beyond where the quasi-isotropic scattered radiation field 
arising, for example, from disk radiation scattered by broad emission-line clouds, dominates the 
point-source disk radiation field is = 0.04i? pc r H) 2 (r sc /0.01)~ 1/2 pc (Dermer and Schlickeiser 1994). 
The appearance of weaker emission line fields in BL Lac objects suggests that accretion-disk rather 
than scattered soft-photon radiation could make a larger relative contribution to the SSC component 
in the inner jets of BL Lac objects, and this component should be included in detailed spectral 
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models using more realistic accretion disk models. The radii where Thomson losses from the 
external scattered radiation field dominates the fields described in the NF and FF regimes are 
derived in Appendices A4 and A5, respectively. Appendix A6 describes different spectral states 
of blazars according to the dominant electron energy-loss processes. Appendix A7 considers the 
transition radius where the electron energy-loss rates from the CMBR field dominates those of the 
external accretion-disk radiation fields. 



9.3. Thomson-Scattered Near-Field Accretion Disk Radiation Spectrum 

Electrons lose energy when scattering soft photons of the accretion-disk radiation field. Con- 
sidering an optically-thick Shakura and Sunyaev accretion disk radiation spectrum, the electron 
energy-loss rate in the Thomson limit of the near field (r ~ R Wr g ) regime is given through 

~ w = 3 caT ( 7 )]r 7 = kNF p- < 41 > 

The scattered vF v radiation spectrum for the near-field component of the disk field is 



f NF^l s 6_^ ( LdJ^ ) ^ N ( ^ /(l+^) g , , 



The derivation of these results is found in Appendix B. 



9.4. Thomson-Scattered Far-Field Accretion Disk Radiation Spectrum 

Again considering an optically-thick Shakura and Sunyaev accretion disk radiation spectrum, 
the electron energy-loss rate in the Thomson limit of the far-field regime is given through 

-IFF = o ccr T h 2~T 2T] r2h 7 a \2 = kpF ^2 ' ( 43 ) 

The scattered vF v radiation spectrum for the far-field component of the disk field is 

- s6 m <l - " 2) ^ ™ ■ ' a ? Jirmfc) ■ (44) 

The derivation of these results is found in Appendix C. By comparing equations (41) and (43), we 
see that the transition from the dominance of the NF to the FF (point source) behavior occurs at 
the transition altitude r< r « 3r 4 r 3 (Dermer and Schlickeiser (1993); see Appendix A3). 
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9.5. Synchrotron Self-Compton Radiation 

The energy-loss rate of electrons as they Compton scatter their self-synchrotron emission is, 
in the Thomson limit, 

-1SSCT =4 -^W 7 u'(e')de'} 7 2 (45) 



3 m e c 2 jq 

The spectral energy density can be directly related to observables through 

/ // /\ rsj /e j cSt var 

e u (e ) = — — , e = — — - and r b ~ — — , (46) 
r£c o 4 (1 + z) 1 + z 

where the last expression relates the comoving size scale to the variability time scale and redshift 
through 5 (Tavecchio et al. 1998). 

The vF v SSC spectrum in the 5-function approximation for the synchrotron and Thomson 
emission spectra is 

<\sA r Jl„.„ ,/ /•min(e',l/e') ,// J 



fe ~ 



(-) 3/2 / de" e"-' N e [(-) l l 2 ]N e [(^) l/2 \ (47) 



(Dermer, Sturner, and Schlickeiser 1997), where e' is related to e in equation (46). Note that a 
factor m e c 2 is missing in the calculation of the SSC spectrum in the paper by Dermer, Sturner, 
and Schlickeiser (1997), so that its importance is underestimated. Tavecchio et al. (1998) give 
expressions when the synchrotron spectrum is approximated by a broken power law. 



10. Model for Blazar Variability 

We operate in the regime where Thomson and synchrotron losses dominate SSC losses, adi- 
abatic losses can be neglected, B can be treated as constant, and processes involving external 
fields can be treated for the thermal accretion-disk model as considered in Section 9. The ana- 
lytic solution presented below can also be extended to the cases where B 2 oc r~ 2 and B 2 oc r~ 3 , 
though we only consdier a constant magnetic field in the calculations presented here. The neglect 
of Klein-Nishina effects on the electron energy loss rate and spectrum is a severe limitation of this 
model (see App. D). Accurate calculations employing realistic spectral forms for synchrotron and 
Compton processes are given in the papers by Bottcher, Mause, and Schlickeiser (1997), Bottcher 
(1999), Mukherjee et al. (1999), and Hartman et al. (2001). 

In the stated approximations, 

-7 = \{ksyn + k T ) + k -^ + ^] 7 2 • (48) 

Location is specified in dimensionless units of f = r/r g , dt' = r g df /(P^Tc), and the coefficients 
follow from the results of the previous section. Solving equation (48) gives 

7 = 7 ( 7i , r, fi) = [7- 1 + F(f, fi)}' 1 , (49) 
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where 111 1 

F(r, n) = [(ksyn + fcr)(r - h) + k FF (-=- - -=) + - ^)] , (50) 

where r is the location at which the particle injection begins. 

Within the framework of first-order Fermi acceleration theory, the injection spectrum of par- 
ticles downstream of the shock is approximated by 

= K^H[ 7l ; lmin , lmax ] , and K = ^Z^H^ (5 l) 

where H[x;a,b] is a Heaviside function such that H = 1 for a < x < b and H = otherwise. 
The term ^f ma x is determined by size scale, radiation, acceleration, available-time and kinematic 
limits (de Jager et al. 1996; Vietri 1998; Rachen and Meszaros 1998; Dermer and Humi 2001). 
The normalization is easily made to the injection power into nonthermal electrons L' e = dE'Jdt' = 
dE*/dt*, where the invariance of injection power between the stationary frame of the black-hole jet 
(starred) system and the comoving (primed) system is invoked. Hence the instantaneous electron 
spectrum at location f is 

Kr /•min[fi,f(t)] 

N>(r,r) = — / dr' Ji~ p H[{~/ - F(r> ,f i )}- 1 ; lmin , lmax ] . (52) 

prlc7 z Jf 

The term n is the location where particle injection ends, and r(t) is the location of the jet that 
emits radiation observed at time t. 

Equation (52) is easily solved using equations (49) and (50) to give the instantaneous electron 
spectrum in the approximation that T, B, and the injection power L' e are constant with time. 
The allowed range of 7 is set by the Heaviside function and the limits on the integral, implying 
T > bmL + ^X^' ^o)] _1 for all values of f > fo- When fo < fi, 7 < [7mi a; + -^ 1 (^ ) ^o)] _1 5 whereas when 
f > r i) 7 < [imax + F{f : fi)]^ 1 . Miicke & Pohl (2000) were the first to modify the the approach 
of Dermer and Schlickeiser (1993) to extended injection. The methods of blast-wave physics (see, 
e.g., Meszaros 2002) may be used to set j m i n and B and treat the spatial evolution of V due to 
interactions with an external medium, at least in the adiabatic and fully radiative regimes. In this 
paper, we simply assign these values. 

Fig. 2 shows the evolution of the spectral components in the toy model for a standard parameter 
set: r = 20, 9 = 1/T, M 9 = 0.1, L Edd = 1, r sc = 0.01, L e = 10 44 ergs s _1 , t var = 1 day, z = 1 
{d,L = 2.4 x 10 28 cm for a cosmology with 70% dark energy and a Hubble constant of 65 km s _1 
Mpc -1 ), B = 1 Gauss, ^ m %n = 10 3 , Jmax = 10 5 , r sc = 0.1 pc, and p = 2.3. Because the electron 
power is small in comparison with the Eddington luminosity, the system could be accreting well 
below the Eddington limit, as argued in an evolutionary scenario by Bottcher and Dermer (2002). 
Different disk models appropriate to smaller values of the Eddington ratio can be treated in future 
work. 



The injection is uniform between 1000 r g and 1500 r g , 2000 r g , 4000 r g , 7000 r g , and 10000 r g 
in Figs. 2a, 2b, 2c, 2d, and 2e, respectively. The relation dr = Pr^c5dt/(l + z) implies 



^3 = 



r(t) 



= 0.8 




— (—? (-) 



(53) 



1000r 9 



where is the observing time in units of 10 4 s, measured from the beginning of the flare when the 
emitting plasma was at tq. For our standard values with Mg = 0.1 and z = 1, a 10 4 s duration 
flare corresponds to the time during which the emitting blob travelling 4000r ff . 

This characteristic time scale is optimum for examining flares with GLAST in terms of flux 
levels for detecting intraday variability and proposed GLAST slewing strategies, as discussed in the 
next section. As can be seen in Fig. 2,, a dominant spectral feature from the near-field accretion 
disk component appears early in the flare, during which the X/j continuum is at a low lever. As the 
flare progresses, the near-field component becomes increasingly weak, and the bulk of the gamma- 
ray emission begins to originate from external scattered soft disk photons. The « 10 MeV - 1 GeV, 
X-ray and synchrotron radiations decline while the > 1 GeV radiation monotonically increases. 
In addition to KN effects on the particle and photon spectra, the diffuse intergalactic infrared 
radiation field will attenuate S> 10 GeV radiation. The important cosmological implications of this 
effect are discussed by, e.g., Primack et al. (1999) and Salamon & Stecker (1998). 

Fig. 3 shows the evolution of the SED with distance from the black-hole core. For constant 
Lorentz factor r = 5 = 20 and Mg = 0.1, the duration of the episode from 10 3 to 10 4 r g is 24 
ksec. The pivoting behavior in the GLAST energy range is evident. Different injection profiles can 
change the the detailed behavior, but the discovery of 7-ray spectral components that individually 
vary in the manner described here would provide evidence for the interaction between the jet and 
the accretion-disk radiation field. At the earliest times in the flare, the actual disk field is apparent 
in the model spectrum at UV, EUV, and soft X-ray energies, and could be revealed from blazars 
at UV/soft X-ray energies during low-intensity states. The direct disk radiation from 3C 279 is 
argued to be detected when 3C 279 was faint in the UV (Pian et al. 1999). 

The EC spectral component associated with the isotropic radiation field would apparently be 
absent or very weak in lineless or weakly-lined BL Lac objects, according to general understanding 
of these sources. In the scenario advanced by Bottcher and Dermer (2002) and Cavaliere & D'Elia 
(2002), BL Lac objects are AGN jet source at a stage in their life where the fueling is in decline 
and the black hole engine is most massive. When GLAST measures BL Lac blazars with better 
sensitivity, limits to the 7-ray flux of a disk radiation component in the 100 MeV - GeV range 
can be used to infer a relationship between the location of the acceleration and radiation sites in 
the jet in terms of S, which can itself be inferred from correlated X-ray and TeV 7-ray variability 
(Catanese et al. 1997) under the assumption that the X-rays and TeV 7 rays have a dominant 
origin in synchrotron and SSC processes, respectively. The near-field and far-field components 
have different antennae patterns, with the far-field component suppressed when observing very 
close along the jet axis (Dermer et al. 1992; Dermer and Schlickeiser 1993). 
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Fig. 4 displays the results of the considered model as a flux density in units of /uJanskys. The 
rapid spectral evolution in the submillimeter is apparent, and can be tested with submillimeter 
and SIRTF observations of the flaring behavior of blazars. The radio spectrum is not modeled, 
and does not properly describe spectra at frequencies less than the mean synchrotron frequency 
( u syn) = 2 x 10 6 B5^ in /(l + z) Hz radiated by electrons with 7 = 7 m j n . The model moreover 
breaks down below the synchrotron self-absorption frequency, which is not treated here. Proper 
modeling of radio requires a consistent treatment of adiabatic losses, which is also not done here. 
Bottcher (1999) points out that the radio spectrum varies differently in flat spectrum radio sources 
and BL Lac objects due to the dominance of the SSC processes in the latter class of sources. 

If a 7-ray flare is a consequence of a relativistic plasma ejection event whereby the plasma 
becomes energized, either through an external or internal shock process as supposed in the model 
approach adopted here, then the delayed emergence of a radio-emitting blob is expected. Evidence 
for this behavior in EGRET data is presented by Jorstad et al. (2001). The quasi-continuous 
monitoring of 7-ray blazars with GLAST, in association with VLBI/SIM searches for superluminal 
motion in radio jet sources, will reveal whether 7-ray flares more typically precede or follow the 
emergence of radio blobs. 

Fig. 5 shows the dependence on observing angle 9 of the SEDs, for the continuous injection 
model where electrons are steadily injected during the period that the blob travels from 10 3 to 
10 4 r s (model le). Note the increasing dominance of the Thomson component over the synchrotron 
component at smaller observing angles 9 due to the different beaming patterns of the two processes 
(Dermer 1995; Dermer, Sturner, and Schlickeiser 1997; Georganopoulos, Kirk, & Mastichiadis 2001). 
Note also that the period of flaring activity for sources observed nearly along the jet is smaller than 
the interval over which the corresponding flux enhancements are observed at large angles to the jet 
axis. These effects are important in statistical treatments of flux- limited samples of blazar observed 
at different angles to the observer direction, which are assumed to be randomly oriented. 

11. Blazar Flare Detectability with GLAST 

We compare our model with the expected sensitivity of GLAST, noting that a recent estimate 
based on the phase 1 EGRET all sky-survey (Fichtel et al. 1994) shows that the rate at which 
GLAST will detect blazar flares sufficiently bright to detect 3a variations on 1 hour time scales is 
about once per month (Dermer and Dingus 2002). 

11.1. Signal and Background in GLAST 

The significance of blazar flare detection with GLAST is estimated (Thompson 1986; Dermer 
and Dingus 2002). The number of source photons with energies E\ < E < Ei detected per unit 
observer time within the solid angle element AO centered around a point source in the direction 
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(9(t), 4>{t)) with respect to the normal of the face of the Large Area Telesccope (LAT) tower arrays 
at time t is 

r E 2 _aq 
S -j E dEA[E,9(tU(t)]<Ps(E,t){l-exp[-^-^]}, ( 54 ) 

where ip s (E,t) is the source photon flux (ph cm~ 2 s _1 E _1 ), and A(E,9,<f>) is the energy and 
angle-dependent effective area of GLAST. We assume that the scattered photons are distributed as 
a Gaussian with A£l u (E) = jTt9^(E) (Thompson 1986), where 9 U (E) is the single photon angular 
resolution. This assumption can be checked against laboratory results to determine the angular 
response of the detector and amplitude of non-Gaussian wings in the point-spread function. The 
GL A S T requirement for the single photon angular resolution (68% containment for on-axis sources) 
is < 3.5° at 100 MeV and < 0.15° at 10 GeV {GLAST Science Requirements Document, 2000). 
Hence 9 U = 0.06{E/E 100 )- 2 / 3 , noting the inverse 2/3 power of the energy dependence of the point 
spread function (Thompson 1986; Fichtel and Trombka 1997). 2 

The source flux ip s (E,t) is related to the vF v flux /^(ergs s _1 ) according to the relation 
Eip s (E,t) = /e/E. We characterize the vF v flux by a power law referred to 100 MeV (= £100) 
photon energy. Thus 

f E = 10- 10 /-io(-^r" ergs cm" 2 s" 1 , (55) 
where a v is the vF v spectral index. 

The energy- and angle-dependent effective area of the GLAST LAT is approximated by the 
function 

A(E, 9) = Aou(9)u(9)(E/E 100 ) a ^ , (56) 

where u(9 = 0°) = 1, a{9 = 0°) = ao, and azimuthal symmetry and time-independence of the 
detector effective area is assumed. An effective area derived from the successful GLAST LAT 
proposal is A w 6200 cm 2 and a « 0.16 for 100 MeV < E < 10 GeV. For 20 MeV < E < 100 
MeV, ao ~ 0.37. This satisfies the requirement for on-axis peak effective area of 8000 cm 2 in the 
1-10 GeV range (GLAST Science Requirements Document, 2001). Hence the on-axis effective area 
A = 6200{E / E 100 ) a cm 2 , where a = 0.4, E < E wo , and a = 0.16, E > E 100 . 

Consider a source whose direction £l s is precisely known. The number of background photons 
with energies E between E\ and E<i detected per unit observer time within solid element AQ of 
the source direction is given by 

B^Afl dE A[E,9(t),^{t)]^> B {E,n s ) . (57) 

The background flux per steradian is denoted by the term &b(E, Q) (ph cm~ 2 s _1 sr _1 E~ v ) and 
is assumed to be time- independent. Considering only high-latitude sources where the extragalactic 



2 The dependence is steeper when E < 100 MeV. 
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diffuse background radiation dominates all other sources of background radiation, 

<S> B (E,n) = K B (-^r a * (58) 
^100 

(Sreekumar et al. 1998) in the range 70 MeV < E < 10 GeV, where K B = 1.72(±0.08) x 10~ 7 ph 
(cnr^-s-sr-MeV)" 1 , and a B = 2.10 ± 0.03. 

An exposure factor X is introduced that crudely takes into account source occultation and 
variation in detector effective area and photon localization with changing source direction of GLAST 
in its nominal slewing mode. For flaring behaviors detected on time scales less than the orbital 
time scale of 90 minutes (£4 = 0.54), X may approach 1. On longer time scales, X ^ 0.2. Hence 

B AnAtXA K B (-^-) 1+a - aB • 100 MeV - 12 Xt 4 (-^) a - 2A3 , (59) 
£-100 1 — 0.4a Eioq 

where we let the solid angle acceptance for background Af2 = AQ U . The approximation in equation 

(59) takes into account that most of the background photons are collected at the lowest energy of 

the range of photon energies. 

Similar approximations and simplications for the number of source counts [eq.(54)] gives the 
result 

The significance to detect a signal at the na level for precisely known background is given by 

n ^7B (61> 

(Li and Ma 1983), where S is the number of source counts and B is the number of background 
counts. 

Equations (59), (60), and (61) characterize high-latitude blazar detectability with GLAST, 
recognizing also that detection requires S > a few counts. At energies E > 100 MeV, a = 0.16, 
and we see that 

£ - t4 ( ^ r ~ a71 ■ (62> 

Detection of blazar sources is favored at « 100 MeV photon energies, except for the hardest sources 
with a u > 0.7 (the mean flat spectrum radio quasar photon index in the 100 MeV - 5 GeV range 
observed with EGRET is (a v ) = 0.2(±0.2) (Mukherjee et al. 1997)). The decline in effective area 
and the increased background at lower energies start to hamper source detection efficiency with 
GLAST when E <C 100 MeV. 



11.2. Light Curves 



Light curves of the model blazar flare shown in Figs. 2-4 are plotted in Fig. 6, except that 
here the nonthermal electron injection ends at r = 5000r 9 . Aftter the radiating blob passes this 
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location, there is no further injection and the nonthermal electrons cool only through radiative 
losses. Thus all three limits specified in the paragraph following equation (52) are relevant. The 
absence of adiabatic losses (see Sikora et al. (2001)), not to mention the assumptions about the 
constancy of B, T, etc., means that this result illustrates only a single limit of parameter space. 
The chief difference in our model from the recent work by Sikora et al. (2001) is the inclusion of 
the direct disk-jet component. 

The high-energy light curves at E = 50 MeV, 500 MeV, and 5 GeV are shown in Fig. 6a. 
The most notable feature is that the 5 GeV light curve continues to harden while the 50 MeV 
and 500 MeV light curves monotonically soften. This is a signature of a disk component fading as 
the blob travels outward, while the component from the scattered disk radiation persists at higher 
energies. The lower energy light curves in this figure reach a plateau after injection stops and electon 
cooling causes scattered accretion disk radiation to dominate in this waveband. The photon energy 
where the two components contribute equally depends, however, on disk parameters such as mass 
accretion rate and inner disk radius, and so might change for different model parameters. 

Fig. 6b shows near infrared, X-ray, and medium energy 7-ray light curves. For this model, the 
synchrotron and SSC components are equal near 1 keV, and the synchrotron component dominates 
at lower energies. Consequently the 0.5 eV IR light curve abruptly declines after the acceleration 
episode ends. The X-ray flaring behaviors seen at 0.5 keV and 5 keV reflect the combined syn- 
chrotron and SSC decline after the injection stops. At 50 keV, the direct accretion-disk component 
causes a late time flattening in the light curve, similar to the behavior seen in the 50 MeV light 
curve, though the latter is due to the scattered accretion-disk component. Note also that these 
light curves will be smeared out if an integration over light travel time effects is made (Chiaberge 
and Ghisellini 1999). 

By comparing the blazar flare light curves at 7-ray energies with the expression for the sig- 
nificance of detection with GLAST, given by equation (62), one sees that the flaring behavior will 
be detectable, though spectral detail will be difficult to extract for the flux given by this model. 
The characteristic signature of a direct accretion disk component will be indicated in GLAST data 
by a soft-to-hard behavior at ~ GeV energies compared to a hard-to-soft behavior at ~ 100 MeV 
energies. 

12. Summary and Conclusions 

In response to a query from Prof. Dr. John G. Kirk about the applicability of the optically 
thin approach used in DS93 to an optically-thick Shakura-Sunyaev accretion-disk spectrum, we 
have reexamined the basis and clarified the approach used in the paper by Dermer and Schlickeiser 
(1993) (a different approach is used in the paper by Dermer, Sturner, and Schlickeiser (1997)). 
An optically-thin formulation is used in the 1993 paper, which appears justified in the case of 
accretion-disk models with large radial flow speeds. For specificity, a photon energy corresponding 
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to the blackbody value was used, though a wide range of choices may be treated. We conclude that 
the method is adequate to treat this special case. 

Energy loss rates and spectral functions were derived in the Thomson regime for the cases 
of soft photons originating from an external isotropic monochromatic radiation field, and from a 
geometrically-thin accretion disk divided into near-field and far-field regimes in the limits r> 1 and 
r 3> Tr g (Appendix A). SSC energy-losses were assumed to be small, and the SSC component was 
calculated in the Thomson regime. A simplified model was examined, and a characteristic variability 
pattern was identified whereby the near-field accretion disk component, initially bright at < 1 GeV 
energies, declines in intensity while the component formed by jet electrons that scatter photons 
of the external isotropic field becomes increasingly dominant at > 1 GeV energies until particle 
injection stops. Campaigns organized around GLAST observations will be crucially important to 
interpret the jet-disk interaction by revealing the physical state of the jet with respect to the disk. 
The intensity of the thermal disk during the low state can additionally be inferred from optical, 
UV, and X-ray measurements of blazars during low-intensity states. 

Both SSC and EC processes contribute to blazar 7-ray production, with a greater contri- 
bution of the internal power dissipated by external Compton (EC) emission in flat radio-spectra 
quasars, and a dominant SSC component in BL Lac objects. Multiple radiation components are 
required to model well-measured contemporaneous blazar spectra (Bottcher 1999; Hartman et al. 
2001; Mukherjee et al. 1999; Sikora et al. 2001). Spectral variations among different radiation com- 
ponents during different intensity states has been considered by Bottcher (2000). A characteristic 
optical/UV spectrum with a = 3/2 is found for injection models with dominating SSC energy losses 
(Chiang and Bottcher 2002). The present paper extends this work to provide an approach to model 
temporal variations of state transitions, though the treatment of SSC losses must be improved in 
future work. 

Our model is very preliminary, and we have not even followed a complete flaring cycle. It is 
unlikely that the magnetic field remains constant over thousands or tens of thousands of r g , and 
analytic models with B 2 oc r~ 2 and B 2 oc r~ 3 , or numerical models with more general variations of 
B remain to be considered, including the related adiabatic losses. The EC component from the soft 
scattered accretion-disk photons will also decline in intensity as the relativistic ejecta leaves the 
broad-line region. Adiabatic expansion will degrade particle energy until a dominant synchrotron 
component remains. The ejecta will expand, cool, decelerate, and on occasion be reenergized by 
internal or external shock interactions. The transition radius from the inner jet, where the external 
radiation field is dominated by either direct or scattered accretion-disk emissions, to the outer jet, 
where the external photon field is CMBR dominated, occurs on the kiloparsec scale for relativistic 
ejecta. 

In the extended jet, the dominant radiation processes are nonthermal synchrotron radiation, 
SSC emission, and Compton-scattered CMBR (e.g., Harris and Krawczynski 2002). Klein-Nishina 
effects can imprint the spectrum with a hardening in the Chandra range, as observed in the knots 



-25- 



of 3C 273 (Marshall et al. 2001; Sambruna et al. 2001; Dermer and Atoyan 2002). As the jet slows 
to nonrelativistic speeds at a terminal shock, either synchrotron or SSC process may dominate, as 
in the western hot spots of Pictor A (Wilson, Young, and Shopbell 2001). Extended jet formation 
could be catalyzed by neutral beam production in the inner jet (Atoyan &; Dermer 2001, 2002). 
The same external-field transformations can be used to solve the problem of the radiation spectrum 
from the outer jet, as presented here for application to the inner blazar jet. 

This paper provides a framework to treat many such problems involving different models and 
geometries of the accretion-disk radiation field in both the optically-thick and optically-thin limits, 
and for other external radiation fields. A simple extension of these results to incorporate advective 
effects on accretion is to change the inner boundary of the Shakura-Sunyaev disk to ~ W0r g , 
within which a low-luminosity or nonradiating advection-dominated flow is found, with convection 
instabilities forming the jet. Application of these results to more complicated accretion-disk models, 
to photomeson production in blazars, and emission spectra in GRBs will be presented in subsequent 
work. 

We thank John Kirk for his interest in our work, and Markus Bottcher for discussions. The 
anonymous referee is thanked for a very constructive report. The work of CD is supported by 
the Office of Naval Research and the NASA Astrophysics Theory Program (DPR S-13756G). RS 
acknowledges partial support by the Bundesministerium fur Bildung und Forschung through DESY, 
grant 05AG9PCA. 



A. Energy Loss Rates in the Near-Field and Far-Field Regimes 

Consider equation (32), 

3GMrh , , (// + Pr f SGMrh 

= 8^f>3 ]_ pT (l-^)3/2 (1+M 4 = 8^5r {lNF + lFF) ' (A1) 

valid in the case [3i = 0, implying an efficiency of 25% for this metric. We have divided the comoving 
energy density into a NF component, with — (5y < fJ-' < 0, and a FF component with < // < ^' max - 

A.l. Near-Field Integral 

The expansion /[/ = —1 + (iV 2 /2r 2 ) in the near-field integral gives, in the limit r> 1, 

f "°jC V (T^fe rJl ^ (A2) 

where 



roo (N 2 -l) 3 
= / dN N 2 (N 2 + l) 4 ~ °"° 23 ' (A3) 
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This result agrees with the estimate q = 2 x 3 3 /(2 2 • 5 4 ) = 0.0216, using the value of the integrand 
at N = T8 = 2 (the integrand peaks at 8' = 1.9/T in the limit T » 1). Fig. 7 presents the integrand 
of equation (A3). 



A. 2. Far-Field Integral 

The far-field integral is 



Iff= L d " (iJwii^ 1 (A4) 



/o ^ (l-M^a + zSr/i') 4 ' 
and ^ maa; = 1 - (i^ in /2r 2 ) in the limit r > R max . Hence 



1 - fr/i 

in the limit a = TR m i n /r <C 1. Therefore 



Mmaz = i = 1 ~ 2a + °( a l V ) . ( A5 ) 



/FF « < 1 4f^- = ^ • ( A6 ) 



A.3. Near-Field/Far-Field Transition Radius 

From equations (Al), (A2) and (A6), 

, 3GMrh ,„ o r , . _. 

implying for a Schwarzschild metric with = 0, the NF/FF transition radius 

f NFFF = T 4 , (A8) 

in normalized units r = r/r g . 



A. 4. Near-Field/Scattered Field Transition Radius 

The location r^x where the NF equals the external scattered radiation field is given by the 
condition 

^ x 21** = U ' NF = u' c = T\l + /f/3H = T 2 (l + ■ (A9) 

Recalling equation (35), we find 



-2/3 -2/3 

f NX = r -f^ = cf 3^/ 3 ^ - (0.59-0.85) ^ , (A10) 

9 T~sc Tsc 
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where the coefficient varies in value as Pi ranges from to 1, and v sc — fgf sc is the characteristic 
size of the scattering region, which we identify with the broad-line region. Hence 

f 2/3 

TNX ~ -f^ ■ (AH) 

Tsc 

Equation (All) holds when f^x < r sc for the simplified geometry of the scattering region assumed 
here, requiring that r sc > l/r sc . If r sc = 0.5 pc around a 10 9 Solar mass black hole, then r sc = 10 4 . 
The scattered radiation will dominate the NF radiation beyond « 2000(r sc /0.01)~ 1 / 3 r 9 from the 
black hole. 



A. 5. Far-Field/Scattered Field Transition Radius 

Equating equations (38) and (43) for the Thomson energy-loss rates of jet electrons by the 
quasi-isotropic scattered disk field and the far field limit of the accretion disk field, respectively, 
gives the location beyond which the direct disk field is no longer important. It is given by 

\/3 f sc . . 

Tsc 

in the limit r» 1. Methods of reverberation mapping, normally applied to nearby radio-quiet 
Seyferts (e.g., Peterson 1993), when applied to flat spectrum radio quasars with strong lines, such 
as 3C 273, 3C 279, PKS 0528+134, etc., can be used to infer t sc and the characteristic size r sc 
of the scattering clouds. A simple dependence of r sc on r g is unlikely, and furthermore, r sc may 
depend upon the spectral state and intensity of the accretion disk. 



A. 6. Intensity States of Extragalactic Relativistic Black-Hole Jet Sources 

These considerations imply a variety of intensity states characterized by different relative values 

of 

-2/3 

-n4 ~ rsc ~ T-i— 2 r sc / a 1 o\ 

tnfff = r , r NX « , r FX = r —q^. (Ai3) 

Tsc Tsc 

These relations hold within the framework of a geometrically thin, steady disk that radiates energy 
through viscous dissipation of gravitational energy in the absence of advective effects. Both the NF 
and the FF spectra are dominant at lower 7-ray energies than the external isotropic component. In 
fact, spectral breaks appear at ~ T 2 lower energy in the FF than the external isotropic component. 

Flaring behaviors for flat spectrum radio quasars could involve transitions from the NF to the 
external scattered field, as modeled in Figs. 2-5. Weaker scattered radiation fields might also involve 
transitions from the FF to the external scattered field. In BL Lac objects, measurements of the 
intensity of scattered disk components with GLAST will restrict the distance from the black hole 
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to the jet, given the specific accretion-disk model. Energy losses due to the SSC component must 
then be properly treated to identify the regime where BL Lac models operate. Modeling within the 
context of an ADAF scenario will improve predictive power and the value of this framework for the 
interpretation of multiwaveband observations of relativistic jet sources. This model is consistent 
with an evolutionary scenario where differences between classes of jet sources are a result of declining 
dust and gas (Bottcher and Dermer 2002), and will be tested through statistics of flat spectrum 
quasars and BL Lac objects (Cavaliere & D'Elia 2002). 



A. 7. Transition from the Inner Jet to the Extended Jet 

A change in spectral state will occur when the pitch-angle averaged electron Thomson energy- 
loss rate due to external radiation fields begins to be dominated by the CMBR rather than the 
accretion-disk emission as a consequence of declining intensity of the accretion-disk field with dis- 
tance. This occurs at the radius r ex j given by 

^u CMB (l + zf = X 47rr2 ^ 4r2 , (A14) 

G.XJ 

where ucmb = 4 x 10 -13 ergs cm -3 , and x an amplification factor to be explained below. 

We write the disk luminosity Ld = rime 2 = ^Edd-^Edd = 1-26 x 10 47 ^EddA^9 ergs s _1 . Solving 
gives 

10fi (x^Edd^) 1 / 2 X 1/2 oon (xlEddM,) 1 / 2 

= 128 T 2 {l + Z? k P c - rxi "T 2 " = 320 (r/10) 2 [(l + ,)/2] 2 PC (A15) 

The radius of the extended jet therefore begins on size scales of order kpc, though with wide 
variation depending strongly on T and less strongly on the accretion-disk luminosity. The quantity x 
represents directional amplification due to the relativistic inner jet. For the synchrotron component 
only (the 7-ray component would mostly be scattered by jet electrons in the KN limit), it seems 
possible that x could reach values of 10 or more along the direction of the inner jet. Deceleration 
of the ejecta to nonrelativistic speeds in distant knots and hot spots would make possible the 
importance of the external Compton component from the inner jet, as defined by the quantity rxj- 
Thus the disk/inner jet radiation field could be important in spectral models of hot spots ~ 10 2 
kpc from the central black hole when the accretion-disk or inner-jet power flares to ~ -^Edd over 
long time scales (S> 10 3 yrs). 



B. Scattered Accretion-Disk Radiation Spectrum in the Near-Field Regime 



From the analysis of the integrand of equation (A2) in the NF regime, we see that the NF 
photons originate from 9' = 1.9/T, or from r « R in the stationary frame. From equations (35), 
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(Al) and (A2), the comoving energy density in the NF regime is 

* = 4^ 4 - 3l+ftcrf( 7 )r2 ' (B1) 

and 4 • 3 1+f3i a = 0.28,0.83 when = 0, 1, respectively. Taking the peak contribution from the 
disk at (i = 1/2, we have the following approximation for the comoving energy density: 

u' NF = Hp 5(e' - Te*/2)8(»' + 1) . (B2) 
The comoving emissivity can be derived from 
j'(e' s X) = c jT 'de' j> dfi' J" 'd 7 j> dn e (l-/3cosV)(^K fe (6 , ,^ , K(7,^e) (^^) (B3) 

(Dermer, Sturner, and Schlickeiser 1997), and cos tp = fi e /j,' + (1 - / Ug) 1 / 2 (1 - // 2 ) 1 / 2 cos(^ e - <fi'). 
Equation (B3) allows to examine the evolution of the particle distribution function in the comoving 
frame. An important simplifying assumption we make at this point is to assume isotropy of the 
electron distribution function in the comoving frame, that is, 

n e ( 7 ,n e ) = ^. (B4) 



For the cross section in the Thomson regime, we make the approximation 

= aT s [e ' s _ 7 V(1 - (3 cos - n e ) . (B5) 

Solving gives 



m , n _ ^ 7 s F „ e(w) , w = ] j . (B6) 

The vF v spectrum 

fe NF = t V b e'J(e' s , Q' s ) , 6 = [F(l - Z^)]" 1 , . (B7) 

a L 

Making the transformation to observer frame quantities using the relation 1 + fi' s — > 5(1 + /j, s )/2T, 
valid when r> 1, gives 



from which equations (41) and (42) follow. A value intermediate to 0.28 and 0.83 is assigned in 
equations (41) and (42). 
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C. Scattered Accretion-Disk Radiation Spectrum in the Far-Field Regime 

The point-source approximation for a monochromatic radiation field coming directly from 
behind the direction of jet motion is 

-Me,n) = ^- c S -^6(e-e). (CI) 

Making use of the relations e = Te'(l + /?r /•*')> M = (/■*' + A")/(l + /?r A*')) and (j> = <j)' gives 

nFF(e ' " } " 8^ r2(i + ^) (C2) 

Following the same procedure as in Appendix B gives 

J (e " ^ } = 32vr^ e -r(l + /3 r ) 7ne(7) ' 7 = [ e "(l - ^) ] ' ( ° 3) 

so that 



/f , = {6 ^ ^ , , = r ys±^ , (C4) 

giving equation (44). 

The accuracy of the expression for j'(e' s ,Cl' s ) in equation (C3) can be checked by deriving the 
comoving electron energy-loss rate 

-m e c 2 7 = j> dn' s J™ de' s Vtf (e' s , fi'J , when V{n e (j) = <5( 7 - 7) , (C5) 

from which follows equation (43). 

D. Klein-Nishina Effects in Compton Scattering 

Georganopoulos, Kirk, & Mastichiadis (2001) have derived an expression for the radiation 
spectrum produced by an isotropic distribution of electrons that Compton-scatters external soft 
photons with monochromatic dimensionless energy e* and energy density u* in the stationary frame. 
Slightly modified, their result is 

_ J3 3^ ,1(1+^ r rf7- to^) HhM2] , (D1) 

where 
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This expression employs the head-on approximation for the Compton scattering process (Jones 
1968; Blumenthal and Gould 1970). The comoving electron spectrum in the range 71 < 7 < 72 is 
given by ^(7). Klein-Nishina effects become important when e > 1/[(1 + z)e*]. 

For a power-law electron distribution N(j) = V^k e j^ p H[y; 71, 72], where V£ is the comoving 
blob volume, they also derive an accurate form for the Thomson regime expression of the specific 
spectral power, given by 

I G kmM) = ^^=(^) (D3) 

U72 ° j [ 46,(p + 3)( 7 25) 2 l+p J W (l+p)(3+p)'' 1 ' 

and show that it reduces to the expression 

LxT (e, n) = vx 2^ (i-jd-pva (D5 ) 

nyp + d) e* 

when p > —1, e <C 4<5 2 e*7|, and 72 3> 71. The Thomson scattered vF v spectrum of an external 
isotropic and monochromatic radiation field was derived by Dermer, Sturner, and Schlickeiser (1997) 
and is given by 

f ,T _ ca T u* !_ Q (l + ^) 2+a r _€_a_ Q ( 

h - a + 2 16tt4 keVb{1 + Z) jj, { ej ' (D6) 

where a = (p — l)/2. Equation (D6) recovers equation (D5) in the limit \x — > 1, though the term 
2(p~ 3 )/ 2 = 0.71,0.84, and 1.0 for p = 2,2.5, and 3, respectively, in equation (D6), is replaced by 
2 p_1 /(p + 1) = 0.67,0.81, and 1.0 for p = 2,2.5, and 3, respectively, in equation (D5). The factor 
(1 + [i)( 2+a ^ I n follows from the range of scattered photon energies and angles in the comoving 
frame, and indicates that fi > is the range of validity of the expression. 

The approach of Georganopoulos, Kirk, & Mastichiadis (2001) avoids the need to transform 
scattered radiation spectra from the comoving to observer frames by directly transforming electron 
energy spectra to the stationary frame. This is useful in the case of specified comoving electron 
distributions. When calculating evolving electron spectra, the external radiation fields must be 
transformed to the comoving frame in order to derive the electron cooling rates. Furthermore, if 
the electron distribution is not assumed to be isotropic in the comoving frame, it is also simpler 
to treat the evolving energy and angle-dependent transformations of the electron spectra in the 
comoving frame. Whether one performs the scattering in the comoving frame and then transforms 
the radiation spectrum, or transforms the electron spectrum to the stationary frame and then 
performs the scattering, the final result must of course be the same. 

Georganopoulos, Kirk, & Mastichiadis (2001) showed that the Thomson approximations (D4)- 
(D6) greatly overestimates the scattered radiation spectra compared to spectra calculated using 
equation (Dl). This discrepancy is considerably reduced for electron spectra that evolve in response 
to Thomson losses in the Thomson-scattering approximation, compared with electrons that evolve 
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under self-consistent Compton losses when using equation (Dl). An approximate analytic treatment 
of KN losses and a detailed numerical treatment of joint Thomson and KN energy losses resulting 
from the CMBR field is performed by Dermer and Atoyan (2002), and Bottcher, Mause, and 
Schlickeiser (1997) calculate blazar spectra using the full Compton cross section. 
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Fig. 1. — Integrand of integral on right-hand-side of equation (32). 
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Fig. 2. — Spectral model of blazar variability, including an SSC component, and EC contributions 
from the direct near-field and far-field disk field and the quasi-isotropic scattered radiation field. 
Parameters and identification of the separate spectral components are given in the legend, (a) 
Electrons are uniformly injected with a comoving power of 1.26 x 10 44 ergs s _1 as the jet travels 
from r = 10 3 r g to r = 1500r 9 . (b) Same as in Fig. 2a, but with the electrons uniformly injected 
between r = 10 3 r 9 and r = 2000r 5 .(c) Same as in Fig. 2a, but with electrons uniformly injected 
between r = W 3 r g and r = 4000r 9 . (d) Same as in Fig. 2a, but with electrons uniformly injected 
between r = 10 3 r 9 and r = 7000r 5 .(e) Same as in Fig. 2a, but with electrons uniformly injected 
between r = 10 3 r 9 and r = 10 4 r s . 




Fig. 3. — Multiwavelength variability behavior for a uniform injection event for the model in Fig. 
2, beginning when the jet is at r = 10 3 r 9 , and ending when r = 10 4 r g . Note the decline of the 
direct disk-radiation field component with distance from the black hole. 




Fig. 4. — Flare profile presented as a flux density in units of micro-Janskys, showing the charac- 
teristic submillimeter/infrared behavior for this class of flare. The radio spectrum is not modeled 
here. 
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Fig. 5. — Angle-dependences of received vF v flux toy model blazar flare. The dominant Compton 
component at small angles reflects the different beaming factors of synchrotron, SSC and EC 
processes. 
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Fig. 6. — Light curves associated with a model flare where injection begins at r = 10 3 r s and ends 
at r = 5 x 10 3 r 9 . Except for the time when particle injection stops, all other parameters are the 
same as in Fig. 2. Time is measured from the start of the flare in units of 10 4 s. (a) Light curves 
at 50 MeV, 500 MeV, and 5 GeV. (b) Light curves at 0.5 eV, 0.5 keV, 5 keV, 50 keV, 50 MeV, and 
5 GeV. 




Fig. 7.— The integrand (N 2 - 1) 3 /[N 2 (N 2 + l) 4 ] of equation (A3). 



